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Abstract 

In the last decade the Sznajd Model has been successfully employed in modeling some properties and 
scale features of both proportional and majority elections. We propose a new version of the Sznajd model 
with a generalized bounded confidence rule - a rule that limits the convincing capability of agents and that is 
essential to allow coexistence of opinions in the stationary state. With an appropriate choice of parameters 
it can be reduced to previous models. We solved this new model both in a mean-field approach (for an 
arbitrary number of opinions) and numerically in a Barabasi- Albert network (for three and four opinions), 
studying the transient and the possible stationary states. We built the phase portrait for the special cases of 
three and four opinions, defining the attractors and their basins of attraction. Through this analysis, we were 
able to understand and explain discrepancies between mean-field and simulation results obtained in previous 
works for the usual Sznajd Model with bounded confidence and three opinions. Both the dynamical system 
approach and our generalized bounded confidence rule are quite general and we think it can be useful to the 
understanding of other similar models. 
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I. INTRODUCTION 



The Sznajd model (SM), proposed in 2000 by Sznajd-Weron and Sznajd, is a model that has 
been successfully employed in the reproduction of some properties observed in the dynamics of 
opinion propagation in a closed community llj]. It is a very simple Ising-like model, that always 
leads to a stationary state of consensus, but with a rich transient behavior. Its originality resides in 
the way the state of the sites evolves: two agreeing sites work together in changing their neighbors' 
state, instead of being influenced by the environment like in the voter model [12l,l3D. 

This model has been extensively studied, either with its original set of rules, or in a variety of 
versions in which one or more rules were changed in order to describe or include specific features 
(see for instance ^4]). It was extended to higher dimensional lattices in 2000 |5, 6] and adapted to 
deal with more than two opinions in networks of different topologies in 2002 [|7D. Also, in 2002, 
Stauffer [js] adapted the bounded confidence restriction, first introduced by Defuant [9] and Krause 
[fioll in models with continuous opinions, to the discrete Sznajd Model scenario. His numerical 



simulation results, for the regular sc 



uare lattice with three opinions, were in disagreement with 



what was found in 2004 by Schulze [|ll|]. that solved the same model in a mean-field approach. 



In this work we propose a new model, with a distinct way of introducing the idea of limited 
persuasion among electors in a discrete model, as the SM. With adequate choices of parameters 
our model restores the set of rules of previous works, allowing comparisons. 

We solved the model's master equation numerically in a mean-field approach and made simu- 
lations of the model in a Barabasi-Albert network and in a square lattice, looking at the transient 
behavior, as well as the stationary state. With the aid of dynamical systems techniques, we were 
able to draw a general picture of its behavior in phase space, identifying its fixed points and basins 
of attraction, for three and four dimensions. This analysis allowed us to understand the origin of 
the discrepancies in the results obtained by Stauffer [0] and Schulze [ nj] for three opinions. We 
also studied the transient behavior of the model for three opinions. Comparing the mean-field time 
series with the ones from the model in a BA network, we see that they share many features, which 
does not happen with the square lattice time series. 

This paper is organized as follows: In the next section we briefly review the rules and the 
behavior of the original SM, the new set of rules introduced by Bernardes and co-authors in 2002, 
the version of the model with bounded confidence introduced by Stauffer in 2002, and finally, 
present our own model. In section III we solve the new model in the mean-field approach, in the 
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general case of M opinions; in section IV we present simulation results, the detailed phase-portrait 
in the special cases of three and four opinions and discuss previous results. Finally, in section V, 
we summarize our conclusions. 



n. THE SZNAJD MODEL 

In the original SM [[it], the sites of a chain with periodic boundary conditions represented voters, 
that could only have opinions (states) cr = ±1. If a pair of adjacent sites had the same opinion, they 
would convince their neighbors with probability p=l; however, if they disagreed their divergence 
would be propagated, with the neighbors adopting an opposite opinion. This model always evolves 
to one of two absorbing states: a ferromagnetic state (consensus state), with all voters with the 
same opinion, or an anti-ferromagnetic state, in which every site has an opinion that is different 
from the opinion of its neighbors (only possible if the chain has an even number of sites or the 
periodic boundary is dropped). The transient, however, displays a rich behavior, that called the 
attention of some physicists 161]. 

This model has been extensively studied, either with its original set of rules, or in a variety of 
versions in which one or more rules were changed in order to describe or include specific features, 
as the possibility of more than two opinions, diffusion of agents, restrictions in the convincing 
capability of agents, or different topologies in the network defining the relationship among vot- 
ers (for a review, see for instance |4]). In most of the works that followed the divergence 
propagation rule was abandoned. 

A. The Sznajd model in complex networks 

In IlZI], Bemardes et al. studied a new version of the SM, that was adapted in order to describe 
the evolution of > 2 opinions in voters located in an arbitrary network. This new model was 
employed to simulate proportional elections with M candidates in a Barabasi-Albert network. In 
their version, each site could be in one of Af + 1 states, the extra state standing for undecided 
voters. Some changes were also introduced in the updating rules, the idea being that, at each time 
step the same average number of neighboring sites were convinced, as in the SM. This can be 
accomplished by setting the probability that a site convinces another one to p = where q is 
the degree (number of neighbors) of the convincing site (in the SM, p = 1 always). Also, a 
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different set of rules was devised for the undecided voters, that were not able to propagate their 
lack of opinion but could be convinced by one of its decided neighbors, even if it did not belong 
to a pair. 

More precisely, the model is defined by the following set of rules: Let a{i, t) be the opinion 
of a site i at time t (a{i, t) E {0, 1, ... , M}, where the positive values represent candidates and 
a{i, t) = stands for undecided voters). Initially, all voters are undecided, except for a set of 
M initial electors (one for each candidate), chosen at random. 

The dynamics consists in visiting each voter in a random (non-sequential) order, applying the 
following rules: 

(I) A voter i is chosen at random. If it is not undecided (cr(z, t) ^ 0), a site j is picked up (at 
random) from the set Fj of neighbors of i and rule II is applied, else nothing happens. 

(Ila) If voter j is undecided (cr(j, t) = 0), then j adopts z's opinion with probability pi = 1/ gj, 
where qi is the degree of site i. 

(lib) If both i and j have the same opinion, voter i tries to convince each one of its neighbors with 
probability Pi = Yg^; 

(He) If i and j have different opinions, nothing happens. 

Like the original SM, this model always evolves towards a consensus absorbing state, but during 
the transient this new model displays a power-law distribution of candidates with v votes. This 
behavior is in agreement with what has been observed in the statistics of real proportional elections 

nnr 
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B. Bounded Confidence 



In all versions of the SM [|15| 1 in which a site i always convinces another one independently of 
its opinion, the system evolves to an absorbing state, in which only one opinion survives. That is 
not always the case in real communities. In an attempt to allow the emergence of different factions, 
Defuant et al. [^, and Hegselmann and Krause [w] introduced the idea of bounded confidence 
for models with continuous opinions (cr g [0, 1]). The idea is to assume that electors i and j, with 
opinions cTj and aj, can interact only if \ai — Cj] < e, that is, if their opinions are close enough. 
This model also evolves to absorbing states, but now eventually with the coexistence of two or 



more opinions that do not interact with each other. This rule can be easily adapted to discrete 
models, as the SM, if opinions are labeled from 1 to M and e is set to 1 . 

In 2002 Stauffer studied the SM with bounded confidence in square lattices [SI], while Schulze 
lUj] studied its mean field version, arriving at different conclusions: In the lattice, Stauffer showed 
that the model 'almost always' evolved to an absorbing state of consensus, while in the mean-field 
approach presented by Schulze the consensus was achieved only in 50% of the cases; the system 
ended up in a state of coexistence of opinions in the other cases. 



C. Generalized bounded confidence rule 



In this work we propose a new model with a generalization of the bounded confidence idea. 
This new model includes, in a single set of rules, the original SM, the complex SM, the SM with 
bounded confidence proposed by Stauffer, as well as many other possibilities. 

In our generalized model, a site with opinion a' has a probability Pa'->a of being convinced by 
another site with opinion a. As in previous models, random non- sequential update is employed, 
and there are no undecided voters. The rules become: 

(I') Choose a voter i at random and a voter j G Fj. If (jj 7^ Uj we do nothing, else we apply rule 

ir. 

(II') Site j tries to convince each one of its neighbors k of opinion aj = ai with probability 
VcTk-^aj ^ coordination of site j. 

Alternatively, rule (11') can also be: 

(II") A neighbor k of j is chosen at random and j convinces k with probability Pa^-^aj ■ 

Note that no assumptions about the probabilities Pa^a' are made beforehand. If Pa^a' is always 
or 1 , then, with a convenient choice of values, one can recover both the usual SM and the discrete 
version with bounded confidence introduced by Stauffer. 



III. TIME EVOLUTION IN A MEAN-FIELD APPROACH 

We can easily write the master equation for the model presented in the last section: 



5 



a' j&k 

- Pa-.a'P{cri = CTj = cr', (Tfc = cr)), (1) 

where Tx and qx are, respectively, the set of neighbors and the coordination (degree) of site X, 
and N is the total number of sites. 
If 

i 

in a mean-field approach, the master equation is reduced to 

cr' 

and in the thermodynamic limit (N oo) we have 

Va = ^{vlr]a'Pa'~.a " vl'VaPa-,a') V CT, (2) 
cr' 

where a time-unit corresponds to a Monte Carlo step (A^ random trials). 

We also define r/ = (771, 772, ... , r]M), F^iff) = t)^ and F = (Fi, F2, . . . , Fm). 
As the sum over a' in (O is antisymmetric with respect to a and a', we get 

f]a = ^ ?7o- = constant. (3) 

cr cr 

From the definition of r] it follows that this constant must be equal to 1 . As a consequence, although 
the flux has M variables, it is restricted to M — 1 dimensions. This also implies that zero is an 
eigenvalue of the Jacobian of F for all values of ff. Also, if ?7cr > V o", the negative term of t),^ 
in (O is proportional to 7]^ (and the term multiplying it does not diverge in the limit r^cr — * 0). So 
the flux is restricted to the region in which all variables are positive (as it should, since 7]^ is the 
probability that a site chosen at random has opinion a in the mean-field). 

The region in phase space where r]^ > and ^7^ = 1 is a regular simplex, and each vertex 
Pi represents the consensus absorbing state of opinion i. The points inside a simplex are unique 
convex combinations of its vertices, suggesting a nice way of representing the phase space of the 
problem. The point P that represents the state (771, 772, 773, . . . , ?7m) is given by 
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p = 

The fixed points of (|2l) are given by 



a 



Va'iVaPa'-^a " Va'Pa~>a') = OT //^ = 0, 

cr' 

and the Jacobian matrix of F, {Jp)a,a' is 

0770-' 



So, in the fixed point, we have: 



Y.a" 'nl"Pa^cT" if ?7a = 

From the expression above it is possible to derive the following conclusions: 

(a) Consider first that a fixed point P* lies in the intersection of manifolds of the type ?7cr = 0, 
so we have (conveniently reordering the variables) 



J* 



V 



(4) 



where is the Jacobian restricted to the non-zero variables in the fixed point and V is 
the Jacobian restricted to the variables equal to zero in the fixed point, which is a diagonal 
matrix. So for each opinion a such that r/* = we have an associated eigenvalue X^r < 0, 



K = -^ivl'fPa^a'- 
a' 

It follows from Q that if x is an eigenvector of J'^ with eigenvalue A, then 
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7^ T 


= A 




V 
















So these eigenvectors and eigenvalues are the same as if we had a model with fewer opinions. 
The eigenvectors are also parallel to all the manifolds r/cr = where P* is. On the other 
hand, if is an eigenvector with eigenvalue \k it follows that J*v^ = \kV^ . So for 
coordinate a we have v^{\k — ^a) = 0. 

Hence, if is not parallel to the manifold defined by r]„ = 0, then we must have A/c = A^r < 
0. 

(b) We focus now on the possible values that A^ may take. If 

r]l,Pa^a' 7^ 0, (5) 

for some cr', then A^ < and the flux, in the neighborhood of P*, is attracted to the manifold 
?7ct = 0. Note that the condition expressed in eq. ([5]) is equivalent to saying that there are 
still sites able to convince a site with opinion a. 

On the other hand, if condition ^ is not satisfied, A^ = and we have, for a point arbitrarily 
close to the fixed point (but outside 77^ = 0), 

a' 

which is the time evolution in second order. Therefore, if there is a a' such that r]*,pa-'^a 7^ 
0, then the manifold ?7cr = is unstable in the neighborhood of the fixed point. If this 
new condition is also not satisfied, then a inevitably interacts only with opinions that don't 
survive in P*, which is the third order term. 

Let O be the set of opinions a such that 77^ = in P*, and let M be the manifold r/^r = 
V a G fi, so P* G M. It follows that, if any opinion in interacts in second order, then 
the trajectories are repelled from M in the neighborhood of P*. If all of them interact in 
first order, then the trajectories are attracted to M. If they all interact only in third order the 
model is degenerate, as opinions in do not interact with opinions outside of it (because of 
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the particular choice of Pa^a')- If all opinions in f2 either interact in first or in third orders 
the model degenerates asymptotically. 

(c) Suppose now that all probabilities Pa^a' are non zero, meaning that all sites have some 
chance of convincing any other one. Consider a surface formed by moving the boundary of 
the simplex inwardly by a sufficiently small amount (but non zero). With the same reasoning 
presented above, we conclude that the flux must come from the simplex's inner region, 
crossing the surface towards the boundary points. 

It follows then that there is an unstable region where all opinions coexist. For M = 3 (3 
opinions) this region is a node. 

(d) Finally, consider a manifold in which all opinions do not interact with each other, that is, 

Va,V<T' 7^0^ Pa'^a = Pa^a' = 0. 

Every point in this manifold is a fixed point, and so, this manifold may have a basin of 
attraction. 

These arguments give a qualitative idea of the evolution of the model in the mean-field ap- 
proach. The less relevant opinions disappear quickly, and the system has a high probability of 
ending up in a state where different and non-interacting opinions coexist (provided these states 
exist). 

IV. THE SPECIAL CASES OF 3 AND 4 OPINIONS 

In the following sub-sections we will analyze in detail the phase portraits that represent the 
dynamics of our model in the cases of three and four opinions, in which they can be drawn. 
We will also present some results about the time evolution of the average number of votes for 
each candidate (the transient behavior) and make comparisons between the mean-field (integrated 
master equation) and the simulated model (in BA networks and square lattices). 
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A. Scenario with three opinions 



It follows from ([3]) that 771 + 772 + 773 = 1, so that the flux is restricted to an equilateral triangle. 
A point in this triangle represents uniquely a set of normalized variables 771,772,^73; the vertices 
Pi, P2 and P3 represent consensus states with opinions 1, 2 and 3 respectively; and the side Aij, 
connecting Pi to Pj, represents the set of states in which opinions i and j coexist. 

One can show that, if Pcr^cr' 7^ V cr 7^ a', (the usual SM lies in this class) the flux has 
an unstable node, where all three opinions coexist; three saddle points, in which two opinions 
coexist; and three stable nodes, representing consensus states (see appendix). Therefore, as far as 
the convincing power among all different opinions is non-zero, the flux is qualitatively the same 
and the model will always evolve to an absorbing state of consensus. However, the basins of 
attraction change, and the same initial condition may belong to different basins if the convincing 
capabilities change (see figure [B. 

But what happens when two opinions a and a' do not interact? This problem, for p^-^cr' = 
j)^,^^ = if |cr — cr'l > 1, has already been studied in detail, both numerically [jSO (square lattice) 
and in a mean-field approximation jlll]. with different conclusions. Stauffer showed that, in a 
square lattice, the stationary state is almost always an absorbing state of consensus in opinion 2. 
However, Schulze simulated the same model in a complete graph - what corresponds to a mean- 
field approach - and only in 50% of the simulations the model evolved to the consensus state with 
opinion 2, found by Stauffer; in the other cases, he observed a steady-state with coexistence of 
opinions 1 and 3. Our analytical approach and generalized model allow us to to understand why. 

In order to understand why the model behaves differently in the mean-field and in the square 
lattice, we first note that the BA network behaves in the same way as the square lattice, almost 
always reaching consensus for opinion 2. One would expect the BA network to behave approx- 
imately like the mean field, as they both have small world properties, unlike the square lattice. 
So whatever process causes the lattice to always reach consensus must also be present in the BA 
network. 

To compare the results for the mean-field and the BA network we integrated numerically the 
equations for the mean-field, to get a phase space portrait of the dynamics. Then, we built an 
'equivalent' portrait for the stochastic model in a BA network, in the following way: we evolved 
the model from an initial condition chosen at random, but with specific expected values of 7/1,7/2 
and ?73, averaging over many simulations. Finally, we plotted the resulting trajectories, together 
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(a) Same convincing strength 




(b) Asymmetric convincing capability 

FIG. 1: (a) Phase portrait for the usual Sznajd Model (everybody convinces everyone with equal 
probabilities), in a mean field approach: there are three stable fixed points (vertices), that corre- 
spond to absorbing states of consensus with opinions 1, 2 and 3; three saddle points, in which 
two opinions coexist, and an unstable node with the coexistence of all three opinions, (b) The 
scenario described in (a) does not change qualitatively as long as the convincing capability be- 
tween any two opinions is different from zero. In this picture pi^2 = P2^3 = P3-+1 = 0.5 and 
P2~>i = P3-^2 = Pi~>3 = 1- The insets resume the interacting rules; the size of the head of an arrow 
indicates the strength of the convincing power in that direction. 
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with the mean-field results (see figure O. 




FIG. 2: Comparison between mean-field trajectories (in gray) and time evolution of the model in 
a BA network (in black), for the case of three opinions and different combinations of Pa^a'- In (a) 
we have e = 1 (usual bounded confidence). The only difference is a slight change in the position 
of the line that separates the two basins of attraction (see detail of central region in (a)). In (b), 
Pi^2 = Pi^3 = P2-»i = P2-»3 = P3-»2 = 1 but ^3^1 = 0, i.e. sitcs with opinion 1 are unable of 
convincing sites with opinion 3. Note that, in both cases, there is no qualitative change, (actually 
little quantitative changes) between the phase portraits of the two networks. 
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The picture shows that in both cases there are basins of attraction for two kinds of solutions: 
consensus in opinion 2 or coexistence of opinions 1 and 3. If the initial opinions are drawn at 
random, with equal probability among opinions 1, 2 and 3, (as done in [ISiiUj]), the initial condition 
will approximately lay in a circle of radius proportional to Yv^, where is the number of sites, 
centered in the point rji = r]2 = rjs = ^/s. 

In the mean-field scenario, this special point is located on the border of the two basins of 
attraction. As a consequence, no matter how large is N and, consequently, how small is the 
neighborhood around the point in which the initial conditions lay, half of its area will be in one 
basin of attraction and half in the other one. 

On the other hand, for the stochastic model this point is, although close to the border, inside the 
consensus basin of attraction. So the coexistence state can only be achieved for small values of N 
(small lattices), for which the fluctuations in the initial condition are bigger. 

For different choices of Pa^a' the qualitative behavior (fixed points, and basins of attraction) 
in phase space is only influenced by which of these probabilities are and which are non-zero. 
When a limit Pa-^a' ^ is taken, typically there will be some fixed points that collapse to already 
existing fixed points where fewer opinions coexist. In the example (b) of figure [2]the saddle point 
between Pi and P3 (coexistence of 2 opinions) collapses to the node in Pi (only 1 opinion) that 
becomes a saddle. 

For all the possibilities where p^^a' is either or 1, the mean-field approach is able to capture 
the whole qualitative behavior of the lattice model (see figure [2(b)] for instance), even though the 
trajectories representing the time evolution of the model in a lattice cross each other, what is 
possible since it is not a flux. This is an assymetric case, for which opinion 1(2) convinces 2(1), 
opinion 2(3) convinces 3(2), but only opinion 3 is able to change opinion 1 (ps^i = 0). 

If we study the time evolution of the average number of votes of each candidate, in the three 
situations studied (square lattice, Barabasi-Albert network and mean-field), we see that the mean- 
filed is a much better approximation for the BA case (see figure [3]) . 



B. Scenario with four opinions 

In the case of four opinions, r/i + 772 + ?73 + ?74 = 1, and the flux is restricted to a tetrahedron. 
With usual bounded confidence rules, Pa^a' = Pa'^a = if |a — cr'| > 1 and Pa~^a' = 1 
otherwise. If we add an interaction between opinions 1 and 4 (pi^4 = p^^i = 1), each one of 
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(c) 



FIG. 3: Comparison between the time series of the Sznajd model in different networks and in 
the absence of bounded confidence. Each color (gray-scale) represents a different opinion. The 
horizontal axis is time and the vertical one is the number (proportion for the mean field) of voters. 
Graph (a) is for a square lattice (approximately 10^ sites), (b) is for a BA network (10^ sites and 
m = 5) and (c) is for the mean-field. We can see the resemblance of the mean field and the BA 
network time series. 

the tetrahedron's faces reproduces the 3 opinion scenario described in the previous section. By 
continuity arguments, we can guess that in this case, there are two distinct basins of attraction, 
shown in figure IH The internal surface isolates completely region I, that includes the edge 
from region II, that includes edge 1 ■<=^ 3. As all points of both these edges are fixed points, there 
are two possible absorbing states with coexistence of opinions (opinions 2 and 4 or opinions 1 
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and 3), regions I and II are therefore the basins of attraction of these states. The fixed point with 
coexistence of four opinions (that lies in the surface between I and II) is unstable, the ones with 
three opinions are saddles (the edges 1<(=^4, 2<(=^3, l-v^2 and 4 3 are unstable manifolds) and 
although the consensus states are in the stable edges (1 <^=^ 3 and 2^4) they are unattainable. 



2 




4 



1 



FIG. 4: (Color online) Boundary between the 2 basins of attraction for the four opinion model 
with e = 1; different basins are in different gray (blue) tones. 



V. CONCLUSIONS 

In summary, we propose a new version of the Sznajd Model, generalizing the bounded confi- 
dence rule. We solve the model in a mean-field approach, for a quite general case, discussing some 
aspects of the dynamics. We showed that the qualitative behavior of trajectories in the mean-field 
approach can be reduced to the study of the cases Pa^a' = and Pa^a' = 1 for each one of the 
possible pairs of opinions a and a'. Also, as long as every opinion interacts with all the others, the 
only possible absorbing state is consensus. 

For the special cases of three and four opinions, that had already been studied in the literature, 
we were able to find a nice way of representing the whole phase space and drew the detailed phase 
portrait, both in a mean-field approach and for a Barabasi-Albert network simulation (in which 
case we developed a method to draw the stochastic trajectories). In both cases the results are 
qualitatively the same (in fact, they are remarkably alike), with two distinct basins of attractions: 
one for an absorbing state of consensus in opinion 2, and another for an absorbing state with 
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coexistence of opinions 1 and 3. The only difference was in the position of the line that separates 
the two basins of attraction. This picture enabled us to understand why in 111 ill (mean-field), the 
model ended up in an absorbing, state of consensus only in 50% of the cases, while in the numerical 
simulations on a BA network IS\\ it almost always ended up in this state. 

Also, regarding the whole time evolution for the average number of electors, with opinion a, 
we were able to derive the following conclusions: (a) for three opinions, a mean-field approach is 
able to reproduce all the main properties of the model when a complex network (usually a network 
with small world properties) is employed to describe the relationship among the electors; however, 
the same does not happen when the model is simulated on a square lattice, (b) The existence of 
any restriction in the convincing power of agents, with at least two opinions that do not interact 
one with the other may lead both in a mean-field approach and in the BA network simulation to 
two classes of absorbing states: consensus in one of the opinions that interacts with all the others 
or coexistence of two (or more) opinions that do not interact. The only difference between the 2 
networks is in the position of the basins of attraction, which means that the initial configuration 
is very important to define the asymptotic behavior, and in order to understand such models, the 
whole phase space must be taken into account. In particular, the 'natural' initial condition with 
a uniform distribution of opinions among the voters may lay in different basins of attraction, if 
different networks (or mean field) are employed. 

The new generalized model introduced by us put the original Sznajd model and a variety of 
bounded confidence versions together in a single model, and the dynamical system's approach 
employed in its analysis allowed to actually understand the whole model and to what extent the 
asymmetries in the way each opinion convinces the others can change qualitatively the behavior 
of the system. We believe that such approach, that is quite general, can easily be adapted to unveil 
new features or draw unifying pictures of other similar models. 



[1] K. Sznajd- Weron and J. Sznajd, International lournal of Modem Physics C 11, 1157 (2000). 

[2] R. A. HoUey and T. M. Liggett, Annals of Probability 3, 643 (1975). 

[3] P L. Krapivsky and S. Redner, Physical Review Letters 90, 238701 (2003). 

[4] D. Stauffer, AIP Conference Proceedings 690, 147 (2003). 

[5] D. Stauffer, A. O. Sousa, and S. M. de Oliveira, International Journal of Modern Physics C 11, 1239 



16 



(2000). 

[6] A. T. Bernardes, U. M. S. Costa, A. D. Araujo, and D. Stauffer, International Journal of Modem 

Physics C 12, 159 (2001). 
[7] A. T. Bernardes, D. Stauffer, and J. Kertasz, European Physical Journal B 25, 123 (2002). 
[8] D. Stauffer, International Journal of Modern Physics C 13, 315 (2002). 

[9] G. Deffuant, D. Neau, F. Amblard, and G. Weisbuch, Advances in Complex Systems 3, 87 (2000). 
[10] R. Hegselmann and U. Krause, Journal of Artificial Societies and Social Simulation 5 (2002). 
[11] C. Schulze, International Journal of Modern Physics C 15, 867 (2004). 

[12] R. N. Costa-Filho, M. P Almeida, J. S. Andrade, and J. Moreira, Physical Review E 60, 1067 (1999). 
[13] M. C. Gonzales, A. O. Souza, and H. J. Herrmann, International Journal of Modern Physics C 15, 45 
(2003). 

[14] F. S. Vannucchi, Master's thesis, Universidade de Sao Paulo - Sao Paulo, Brazil (2006). 

[15] Without the divergence propagation rule. In the original version a paramagnetic state was possible. 

VI. APPENDIX 

In the fixed point we have 



From this set of equations its is trivial to show that, if only one of the r]cr 7^ 0, we have a stable 
fixed point in one of the vertices; If i]^ ^ for all values of a, we can define /ii = vi/r^^ and 
= ^^/v3', equations ^ can then be written as 



(6a) 



VhlPl^2 - rilri2P2-,l + V2V3P3^2 - V2ViP2^3 = 0. 



(6b) 



(6c) 



/Ul/U2P2^1 - f4Pl-2 + AtlP3->l - Pl~>3 
/Xl/i2Pl^2 - iAP2^1 + /^2]?3^2 - P2^3 



0. 







(7b) 



(7a) 



from (fTbl) we get 



filp2^1 + P2^3 
yUlPi^2 + P3^2 ' 



17 



and (|7al) becomes 



/ilP2-.l(/i?P2^1 +P2^3)(AHPi^2 " Pl^2(/i?P2->l + ^2^3)^ + 

+ (/iiP3^1 -pi^3)(/iiPi^2 +P3^2)^ = 0, 

that is an ordinary polynomial of third order in fxi: 



/(/il) = {p1^iP3~>2 + P3^iPL2)/^1 - Pl->2(P2-.lP2^3 " '2P3~>lP3^2 + 

+ pi^2Pl-,3)f4 + P3^2{P2^lP2^3 " 2pi^2Pl^3 + P3^lP3^2)/il- (9) 

- (P1-2PL3 + Pi^3p1^2) = M + Bfij + Cfii + D = 0. 

The real positive roots of this polynomial corresponds to fixed points in which all three opinions 
coexist; if Pa-^a' > 0, A > and D < 0, there is 1 or 3 positive roots. 

Suppose, by absurd, that there were three real positive roots: we would then have S < 0, C > 
and the discriminant A = AAC^ + AB^D - B^C^ + 27 A^D^ - 18ABCD < 0; 
defining 

C , . B 

a = r and p = t-, 

3{AD^)^ 3{AW)^ 

we note that, if A, C > and B,D < 0, we have a,(3 > 0. Besides, = 27a^ AD^ , B^ = 

27(3^ A^D and BC = 9ADa(3, 

so that 

A = lOSa^A^D^ + 108p^ A'^D^ - SWp^A^D^ + 27A^D^ - 162a(3A^D^ 

and 

S = l^f^ = + 4/?^ + 1 - 3a'/32 - 6a/5 < 0, 

that is 



5{a, p) = E{p)a^ + F{p)a^ + G{(3)a + H{(3), 

where E = A,F = -S/?^, G = -6/3 and H = Ap^ + 1 (note that H > 0). So, if we fix p, S{a) 
has a negative root. But E > 0, and, unless 5 has a positive root, we would have 5 > for a > 0; 
However, this is only possible if 6 has three real roots, what is equivalent to 
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AEG^ + AF^H - F^G^ + 27E^H^ - 18EFGH < 

or 

- 3456/5^ - 108/3^(4/5=^ + 1) - 324/5^ + 432(16/3^ + 8(3^ + 1) - 1296/3^(4/5^ + 1) = 
= - 432(/5^ - 1)2 < ^ /5 > 1. 

But because 5{a, (3) = 6{(3, a), we have a>l^a(3>l^ BC < 9AD. However, we also 
have 



BC - 9AD = 8pl^^Pi^2P3^2Pl^3 + ^P2~*iPI-^2Pi^3 + QP3^lPl^2P2^3 + 
+ 4:P3_ip1^2P1^2Pi^3 + P2-.iP1^2P3^2P1->3P2~.3 + 

+ P2^iP3^iPi^2p1^2P2^3 + 2pl_,^pi^2pl-,2 + '2pl_,2P3-^2pl^3 > 0' 

what is a contradiction! So, there is one and only one positive root, and only one (unstable) fixed 
point inside the triangle. The corresponding values of t] can be obtained from (l6al) : 

Vi = : — V2 = ; — r and 773 



A^l + /i2 + 1 /il + Ai2 + 1 A^l + yU2 + 1 

Finally, it is also easy to see that, if only one of the 77 are equal to zero in we have three 
other (stable) fixed points: 





= 0, 


^2(3) 


V2 


= 0, 


V3{1) 


V3 


= 0, 





P2(3)^3(2) 


P2(3)- 


-*3(2) + P3{2y 


-2(3) 




P3(l)-1(3) 




P3(l)- 


--1(3) +Pl{3)- 


-3(1) 




Pl(2)^2(l) 




Pl(2)- 


-*2(i) +P2{iy 


-1(2) 



and 



By continuity reasons there must be saddle points within the unstable manifolds along the sides. 
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